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The normal state single particle spectral function of the high temperature superconducting 
cuprates, measured by the angle resolved photoelectron spectroscopy (ARPES), has been considered 
both anomalous and crucial to understand. Here, we report an unprecedented success of the new 
Extremely Correlated Fermi Liquid theory by Shastry to describe both laser and conventional syn- 
chrotron ARPES data (nodal cut at optimal doping) on Bi2Sr2CaCu20s+(5 and synchrotron data on 
Lai.85Sro.i5Cu04. It fits all data sets with the same physical parameter values, satisfies the particle 
sum rule and successfully addresses two widely discussed "kink" anomalies in the dispersion. 

PACS numbers: 71. lO.Ay, 74.25. Jb,74.72.Gh,79.60.-i 



Angle resolved photo-electron spectroscopy (ARPES) 
was the first probe to provide a detailed view of the 
anomalous nature of high temperature cuprate supercon- 
ductors, discovering unexpectedly broad spectra with in- 
tense and asymmetric tails that have remained an endur- 
ing mystery for the last two decades. Conventional data 
taken with high energy (> 15 eV) photons from syn- 
chrotron light sources have recently been supplemented 
with laser ARPES data [l|, l2| from lower energy (6 or 7 
eV) sources. The latter show considerably sharper fea- 
tures near the Fermi energy. A drastic possibility to ac- 
count for this distinction is that the sudden approxima- 
tion could break down for the smaller photon energies 
used in laser ARPES Q 

An important un-answered question is whether the re- 
sults of the two spectroscopies could be reconciled in a 
single theoretical framework that does not abandon the 
sudden approximation. More broadly, can we understand 
the wide variety of observed lines shapes in a theoretical 
framework with a sound microscopic basis and a single 
set of parameters? 

In this Letter, we confront a recent theory of Extremely 
Correlated Fermi Liquids (ECFL) proposed by Shastry 
Q with the above challenge. The new formalism is com- 
plex and requires considerable further effort to yield nu- 
merical results in low dimensions. In the limit of high 
enough dimensions, however, a remarkably simple ex- 
pression for the Green's function emerges; it is signifi- 
cantly different from the standard Fermi Liquid Dyson 
form, while satisfying the usual sum rules. We use this 
simple version of ECFL Green's function in this Let- 
ter, motivated by the attractive spectral shapes produced 
with very few parameters Q • In this Letter we show that 
already the simplest version of the ECFL theory, with 
very few parameters, is successful to an unprecedented 
extent in detailed fitting of a wide variety of normal state 
cuprate ARPES line shapes. Interesting predictions are 
made for the higher temperature spectral line skew. 

Our focus in this Letter is on the data of optimally 



doped Bi 2 Sr 2 CaCu20 8+ 5 (B12212) and Lai.gsSro.isCuCU 
(LSCO) superconductors in the normal state, taken 
with k along the nodal direction connecting (0,0) to 
(7r/a, ir/a). Most of the data is taken from the published 
literature, while some original data are also presented 
(Bi2212 data in Figs. I4|5p . Our sample is an optimally 
doped Bi2212 (T c = 91 K), grown by the floating zone 
method at the Brookhaven National Laboratory (BNL), 
and was measured at the Stanford Synchrotron Radiation 
Lightsource (SSRL) beam line 5-4 using 25 eV photons. 
The resolutions are 15 meV (energy) and 0.3° (angle). 

Line shape model: The ECFL spectral function is 
given as a product of an auxiliary Fermi Liquid (aux- 
FL) spectral function AFL(k,uj) and a second frequency 
dependent "caparison" factor 0, H| : 
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where n is the number of electrons per Cu02 unit cell, 
(X)+ = max(X,0), f g = (l - f ) e(k), where e(k) is 
the bare one-electron band dispersion (see later). Here, 

with 



A FL (k,uj) = 1 -^m^^; 
<£>(cj) = 



uj 2 + r 



$(u;) 
2 

■ exp 



UJ 2 +T 2 



(2) 



where r = 7r/c#T, T is the temperature, and uj is to be 
understood as uj — i0 + . Here, ujo is the aux-FL energy 
scale (i.e. high uj cutoff), and governs the lifetime, 
and, by causality, the quasi-particle weight (i.e. the wave 
function renormalization) of the aux-FL, Zfl = (1 + 
-^f^)~\ as identified from Ke^Q. 

The ECFL energy scale Ao measures the "average in- 
trinsic in-elasticity" of the aux-FL. It is given [4[ as 



dwf(uj)(A FL (k,oj)(^-uj)) 
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where (-)bz denotes averaging over the first Brillouin 
zone. 
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FIG. 1. (a) A as a function of c^o for various Zfl- Other 
primary ECFL parameter values are n — 0.85, T — 100 K, 
and ^ as described in the text. A small n value, 0.010 eV, was 
used for this plot, which is used as a "lookup table" during the 
fit. (b,c) Examples of the spectral function calculated with 
different values of the effective sample quality parameter r\. 
See the caption of the next figure for parameter values used. 
The instrumental energy broadening of 10 meV (FWHM) is 
included. 



The parameters that enter this description are now 
listed. The "primary parameters" defining the ECFL fit 
consist of the dispersion £g taken from band theory, the 
density n, temperature T, and the aux-FL parameters 
Aq, co>o, Zfl-, ^o- Of the last four parameters, only two 
are free parameters. For instance, c^o and Zfl can be 
taken as free parameters, and VLq and Ao can be calcu- 
lated using the equation for Zfl and Eq. [3j respectively. 

The parameter n in Eq. (|2]) is an additional "secondary 
parameter" Q with respect to the ECFL theory Q. Its 
origin is in impurity scattering as argued in [8], and ad- 
ditionally, in scattering with surface imperfections. Our 
fits determine n ~ 0.03 eV for laser ARPES and n ~ 0.15 
eV for conventional ARPES. Greater penetration depth 
in laser ARPES suggests that it should be less sensi- 
tive to surface imperfections, thereby yielding a smaller 
f]. We therefore propose that this parameter summarizes 
the effective sample quality in different experiments. The 
difference in line shapes arising from these values of 77 is 
demonstrated in Figs.[TJb,c). 

Our strategy is to fix a common set of intrinsic param- 
eters for all the materials, and allow rj to be determined 
separately for each class of data. The most time consum- 
ing part is the calculation of Ao, the results of which are 
summarized in Fig. [TJa). 

In our line shape analysis (1) we first set n = 0.85, cor- 
responding to the optimal doping. (2) Here £g is taken 
to be the un-renormalized band dispersion, taken from 
the literature [9] , and then scaled to fit the observed oc- 
cupied band width, 1.5 eV, of the Bi2212 ARPES result 
[10] [11]. (3) We choose Zfl — 1/3, to account for the 
dispersion renormalization due to the high energy kink 
[id . Il2j . which in this theory is caused by the energy 
scale coo (cf. Fig. [5j). (4) Finally, in all simulations, we 
include the finite energy resolution effect and the finite 
angle resolution effect as a combined Gaussian broad- 
ening (10 meV FWHM for laser ARPES and 25 meV 
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FIG. 2. Laser ARPES data (symbols, B12212) from Ref. [3 
fit with the ECFL line shape (red lines) . The free parameters 
of the fit were cjo (0.5 eV), n (0.032 eV), and ^ (shown). 
Fixed parameters: n (0.85), Zfl (1/3). Derived parameters: 
A (0.12 eV), Q (0.14 eV). Other than n and ^, the same 
parameters are used elsewhere in the Letter. In (a), the gray 
line corresponds to the theoretical curve with £g = 0.15 eV. 



FWHM for conventional ARPES) in energy Jig. 

Line shape fit for laser ARPES: Fig. [2] shows the 
fit of the laser ARPES data with the ECFL line shape. 
These fits were made using a procedure that is somewhat 
more restrictive than that in the recent work of Casey 
and Anderson [3, EH , using the equations of Doniach- 
Sunjic, Anderson-Yuval, and Nozieres-de Dominicis [l6| 
(CADS): we are using global, rather than per-spectrum, 
fit parameters. However, our fit is somewhat less re- 
stricted than other fits shown in this Letter: here we 



allow a small variation of £g as in Ref. Il4j. We find an ex- 
cellent fit quality, at least comparable to CADS [14]. The 
gray line in panel (a) shows our calculation for k > kF- 
Our expectation is that, were the data for k > kp avail- 
able, we would find a reasonable fit in this k region as 
well [17], as for other data sets below. 

Line shape fit for conventional ARPES: We find 
that the magnitude of the parameter ujq (0.5 eV) deter- 
mined from the fit of the sharp laser data works very well 
also for the conventional ARPES data [19]. Thus, all pa- 
rameters other than 77 are fixed, with one small exception 
in Fig.[4t[d), where a slight change in ujq produces a much 
better fit over a larger energy range for LSCO. 

Fig. [3] shows our fit of the data in Ref. [l8| with a single 
free parameter n. The amount of the "extrinsic back- 
ground" (bg) in ARPES is an issue of importance [2ol — 
l22[, especially when analyzing the conventional ARPES 
data. Here we fit the bg subtracted data, as well as the 
raw data (panel d). For subtracting the bg, we use an 
often-used procedure [H, Q of equating the background 
to a fraction ("bg scaling factor") of the data far be- 
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FIG. 3. Conventional ARPES data (Bi2212) fit with the 
ECFL line shape. The data are from Ref.[H (T c = 90 K). (a) 
The data (symbols) and the fit (red lines) are shifted verti- 
cally by the same amount for ease of view, (b) An example of 
the raw data and the fit data is shown for ki. The background 
(bg) spectrum (see text) was subtracted from each raw data, 
and the resulting data, shown in (a), are then fit. (c) The 
fixed parameters used for the fit. Thus, in this figure, r\ is 
the only fit parameter (cf. Fig. [2] caption), (d) Raw data at 
k — Uf fit with a somewhat greater 77 value, (e) The current 
fit compared with a fit using the CADS line shape. 
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FIG. 4. Conventional ARPES data, including our own (a,b), 
fit with the ECFL line shape. The procedure used to fit these 
data are identical with those of the previous figure, i.e. a fit 
with a single free parameter 77, with (d) being a single ex- 
ception, (a, b) Optimally doped Bi2212 (T c = 91 K). (c) 
Optimally doped LSCO data [27|,[28|. (d) A test fit up to 0.6 
eV for the LSCO data with a small change to ujq for the same 
data as in (c) but over a wider energy range. By changing 
ujo slightly from 0.50 eV to 0.42 eV, we see that an excellent 
fit up to 0.6 eV is found. The LSCO data, as far as we are 
aware, is fit only by the ECFL theory, since an energy depen- 
dence rising linearly for occupied states occurs naturally and 
uniquely in the ECFL spectral function. 



yond the Fermi surface crossing (k = fcio for this data 
set). The bg scaling factor, 1/2 for this figure, is deter- 
mined to be the maximum value for which the resulting 
intensity is not negative. As shown in the panel d, the 
ECFL fit remains good by adjusting 77, whether or not 
the extrinsic background is subtracted. In contrast, we 
find that the CADS theory, notwithstanding its notable 
successes [14|, UM •> cannot cope with even the background 
subtracted data (Fig. 3e), giving too steep a fall off to- 
wards the left. Likewise, the MFL fits (8ll25| have been 
shown to compare well with the data only after substan- 
tial background subtraction [H, [26| . 

Our own data on Bi2212 data, taken at T c and well 
above T c , covering a similar temperature range as the 
laser data of Fig. [2j can be fit equally well with the same 
background subtraction procedure, i.e. with the "bg scal- 
ing factor" (1/2), as shown in Fig. 3J 

We also find that the data for a lower- T c cuprate LSCO 
can be fit very well with the same intrinsic parameters. 
Here, we shall discuss only the k — kp data for brevity. 
In this case, we determine that the "bg scaling factor" be 
1. The subtracted "bg" data 13, HH is shown as the gray 
curve in Fig. 0Tc). Given their weak superconductivity 



features 27|, [28|, these LSCO data are taken to repre- 
sent the normal state property even if the temperature is 
slightly lower than T c . As for the Bi2212 case, the data 
can be fit well even without the background subtraction, 
if a somewhat greater 77 value (« 0.17 eV) is used. It 



is clear, from Fig. lU^c), that the data at a temperature 
as low as 25 K can be fit very well with the ECFL line 
shape. In addition, in working with LSCO line shapes, 
we noticed a steady and rapid rise in intensity beyond 
— uj = 0.25 eV, a behavior different from that of Bi2212. 
We leave the full discussion of this non-universal behav- 
ior for future work. However, we find it exceptional that 
the current theory is able to describe the line shape of 
LSCO up to very high energy, as shown in Fig. 11(d). 

Kinks in the spectra: The two independent energy 
scales uo and Ao are determined from our fit as ~ 0.5 
eV and ~ 0.1 eV. These are natural candidates for the 
two main dispersion anomalies in the cuprates [H, Hsl as 
in Fig. [5] Well-defined energy distribution curve (EDC: 
intensity curve at a fixed k value) peaks disappear in 
a wide energy from ~ 0.3 eV to ~ 1 eV, as observed 
experimentally for the high energy kink [10|, LL2| • As this 
feature already exists in aux-FL, it cannot be associated 
with Ao but rather with cjo- The (numerical) dynamical 
mean field theory [30| can already account for this feature 
as can the present ECFL (analytical) theory. 

Turning to the low energy ARPES kink at ~ 70 meV, 
Figs. [5fc,d,e) illustrate the observed weak dispersion 
anomaly in the normal state data (c), reproduced in the 
ECFL theory (d) but not in the aux-FL theory (e). Here 
we use a visualization method for momentum distribu- 
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FIG. 5. Image plots of the spectral function for (a) the ECFL 
theory (77 = 0.17 eV), and (b) that of the auxiliary FL. (c) 
The data of Fig. HI a) before (inset) and after (main) MDC 
normalization, by which each MDC is scaled and shifted to 
have minimum (green) and maxmimum 1 (red). Blue cor- 
responds to 1/2. (d) The near- Ef part of the ECFL spectral 
function of (a), after MDC-normalization with MDC peak 
positions traced by black line. The black line in (c) is from 
(d). (e) The near-^F part of the aux-FL spectral function 
of (b) after MDC-normalization. The MDC peak positions 
are traced by gray dashed line, while the black line is from 
(d). The bending such as shown by the black line here is 
commonly referred to as the kink, (f) The temperature de- 
pendence of the peak asymmetry compared for three different 
theories, rising as T 2 for ECFL. Theory parameters for the 
calculation, apart from T, are taken from the fit of Fig. [2{a) 
for the ECFL, from the equivalent fit of Ref. H for the CADS, 
and from Ref. 26 for the MFL. 



tion curve (MDC: intensity curve at a fixed uj value), an 
object discussed primarily for low energy kinks. Thus 
this feature originates from the scale Ao, it causes an in- 
creased asymmetry and the (blue) shift of the peak to 
high hole energy, when the third term in the caparison 
factor (x^a^~) °f E becomes important. To our 
knowledge, the ECFL theory is a unique analytical the- 
ory that has both these kink features arising from purely 
electronic (extreme) correlations. 

In Fig. EJf), we show the temperature dependence of 
the dimensionless peak skew or asymmetry, defined as 
(HL - HR) / (HL + HR), where HR (HL) is the half- width 
at half maximum on the right (left) side of the peak. 



The predicted T-dependent asymmetry, predicted even 
greater for 77 « 0.15 eV (synchrotron data; not shown), 
would be interesting to explore in the future. 

Further work is necessary to refine the picture sug- 
gested in this Letter. For example, as — £g increases, the 
line shape becomes somewhat too asymmetric. Work is 
also in progress to apply the theory to two particle re- 
sponse as seen, e.g., in optical conductivity. We have 
checked that the bubble approximation (conductivity as 
a product of two G's) shows an agreement in the order of 
magnitude of the frequency scale and the conductivity. 

Conclusions: We have shown that it is possible to 
understand both ARPES data sets (laser or conven- 
tional) comprehensively, with identical physical param- 
eters. Work going beyond the nodal cut and the optimal 
doping value is in progress. The theory is very toler- 
ant of the uncertainty in the background subtraction for 
the conventional ARPES data. Additionally, the theory 
satisfies the global particle sum rule, and contains two 
inter-dependent energy scales (ujo and Ao) that corre- 
spond well to the energy scales of the two kinks. Thus 
the simplest version of the ECFL theory using a small 
number of parameters, provides a framework to under- 
stand the ARPES line shape data for the normal state of 
the cuprates: it works extremely well across techniques, 
samples and temperatures. 
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